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I derive the overlap Dirac operator starting from the overlap formalism, discuss the 
numerical hurdles in dealing with this operator and present ways to overcome them. 

1. Introduction 

The overlap Dirac operatorS, derived from the overlap formalism! for the special 
case of vector gauge theories, is a way to realize exact chiral symmetry on the 
lattice. Exact chiral symmetry on the lattice does come at a price - numerical 
implementation of the overlap Dirac operator is significantly more expensive than 
Wilson or staggered operator. In spite of this numerical hurdle, we already have 
several physics results in quenched gauge theories using the overlap Dirac operator: 

(i) Evidence for spontaneous chiral symmetry breaking at zero temperature!. 

(ii) Evidence for chiral symmetry breaking in the deconfined phase possibly due 
to a dilute gas of instanton and anti-instantonsu. 

(iii) Evidence for a diverging chiral condensate in the two dimensional U(l) casefx 

(iv) A study of exact zero modes of overlap fermions in the adjoint representation 
lend some support to the existence of fractional topological chargeQ. 

In this talk, I shall derive the overlap Dirac operator starting from the overlap 
formalism, discuss the numerical hurdles in dealing with this operator and present 
ways to overcome them. 

2. The Overlap formalism 

The determinant of the chiral Dirac Operator C = cr^d^ + iA fl ) can be realized 
on the lattice as an overlap of two many body statesl'i, namely; 

detC=< 0- |0+ >, (1) 

where |0± > are many body ground states of a^H(m)a and a^j^a respectively. 
The a' and a are canonical fermion creation and destruction operators and 75-Zf (to) 
is a massive Dirac operator on the lattice with the mass set to a value less than 
zero. One choice is the Wilson Dirac operator, H(m) = H w (m). This realization of 
the chiral Dirac operator is natural since C is an operator that maps two different 
spaces, namely spinors under the (0,1/2) representation to (1/2,0) representation. 
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Therefore C does not have an eigenvalue problem and the determinant of C is a 
map between the highest form in the two spaces connected by the operator C . 

Clearly, the overlap formula does not fix the phase of |0+) since it is only defined 
as an eigenvector of a Hamiltonian and this is how it should be since the chiral 
determinant is a map between two different spaces. The details involved in the 
phase choice and possible gauge breaking is the subject of chiral gauge theories. 
For vector gauge theories, we want det CC^ — |(0 — |0+)| 2 and the phase choice 
does not matter indicating a trivial cancelation of anomalies. 

Computing the overlap of two many body states seems like a insurmountable 
numerical task in four dimensional theories since one has to diagonalize H w , form 
the many body state from the negative energy single particle states and compute 
the overlap by computing a determinant of a dense matrix, half the size of H w . But 
there is an elegant solution to circumvent these steps by. directly dealing with the 
many body states and this is the overlap Dirac operatora. 

The massless overlap Dirac operator 



D =± 
2 
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is derived from the overlap formalism as follows. Let U be the unitary matrix that 
diagonalizes H w : 



H W U = U\; U=[Z J , . 
Using det U = det a/ det 6^ , we derive 
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It is not immediately clear as to how it helps numerically since one will have 
to deal with e(H w ) without having to diagonalize H w . There are two possible 
approaches. One approach is to use Gegenbauer polynomials to represent y^^pm. 
Typically one need to go to a high order polynomial and this method is Jiot expected 
to efficient. The other approach is to use the rational approximationE2l where one 
approximates e(H w ) as a sum of poles 



e(H w ) = c + 



CjH u 
Hi + 



(5) 
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Using the method of multiple masses, one action of e(H w ) on a vector can be realized 
by a single conjugate gradient algorithm independent of the number of poles. This 
makes it numerically quite attractive. 

3. Spectrum of the quenched H w 

e(H w ) is discontinuous at the zero of H w . Approximations have to be good up to 
lowest eigenvalue of H w and this can be a problem if H w has very small eigenvalues. 
The density of the spectrum of H w (m), p(X), in a quenched ensemble has a non-zero 
p(0) at any fixed lattice coupling at the values of m that are relevant (m < m c )Ll 
This numerical result has support from an analytical argument where one shows 
that small defects can already give rise to a gapless spectrumEx 

One can also show that a change in gauge field topology necessitates zero eigen- 
values at any mass. To see this, let us assume we have a gauge field configuration 
that has zero topology. Then H w (m) has an equal number of + and - eigenvalues. 
Consider evolving from this configuration to another gauge field configuration that 
has non-zero topology. This configuration has a spectrum where the number of + 
and - eigenvalues of H w {m) are not equal. The spectrum as a function of the evo- 
lution has one configuration in the path where H w (m) has a zero eigenvalue. In a 
discrete evolution scheme the exact zero will be avoided but one can have arbitrarily 
small eigenvalues. 

Therefore one will have to live with very small eigenvalues of H w (m) or its 
variants. Numerical techniques that deal with e(H w ) will have to project out a 
few small eigenvectors and treat them exactly. On a finite lattice and at a fixed 
lattice spacing, the number of eigenvalues below a fixed number A m i n will grow with 
volume since p(0) is finite. This would mean that one has to project out more 
eigenvalues as one increases the volume and/or go to a larger number of poles in 
the rational approximation. It is useful to compare the overlap formalism with the 
related method used to realize chiral symmetry on the lattice, namely domain wall 
fermionsE3. This is a five dimensional realization and the effective overlap Dirac 
operator is obtained by setting H = Hd = \og(T w ) where T w is the transfer matrix 
in the fifth directions. The low lying spectrum of is completely governed by 
the low lying spectrum of H w and hence the problems caused by a finite p(0) exist 
for domain wall fermions0. In practice one works with a finite extent in the fifth 
direction (L s ) and this amounts to an approximation of the e(H w ) by t&nh(^L s Hd). 
Clearly, small eigenvalues are not taken care of properly at a finite L s and one will 
have to go to larger L s as one increases the lattice-nplume at a fixed lattice spacing. 
Current simulations using domain wall fermional3 seem to indicate a significant 
effect due to finite L s . One can avoid this by projecting out small eigenvalues and 
treating them exactly in the domain wall formalism!!!!. 

Each action of e(H w ) requires a Conjugate Gradient type algorithm and therefore 
the solution to the equation of the form D (m)ip = b requires nested Conjugate 
Gradient. Is it numerically much more involved than domain wall fermions since it 
only involves one inversion of a higher dimensional operator? One can write down a 



4 Recent Results Using the Overlap Dirac Operator 



five dimensional operator from which one gets the required four dimensional overlap 
Dirac operator by integrating out all but one fermion degree of freedom. An analysis 
of the condition numbers shows that the the five dimensional inversion is no less 
expensive that two nested conjugate gradientsli3. In the nested case, it is easy to see 
that the condition number is proportional to the product of the condition number 
of H w and the fermion mass, /i. This also turns out to be the case for the five 
dimensional version and for the conventional domain wall fermions. This shows 
that it is practical to work directly with the four dimensional operator. 
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